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ABSTRACT 

Galaxy clusters, the most massive collapsed structures, have been routinely used to determine cosmological 
parameters. When using clusters for cosmology, the crucial assumption is that they are relaxed. However, 
subarcminute resolution Sunyaev-Zel'dovich (SZ) effect images compared with high resolution X-ray im- 
ages of some clusters show significant offsets between the two peaks. We have carried out self-consistent 
N-body/hydrodynamical simulations of merging galaxy clusters using FLASH to study these offsets quantita- 
tively. We have found that significant displacements result between the SZ and X-ray peaks for large relative 
velocities for all masses used in our simulations as long as the impact parameters were about 100-250 kpc. Our 
results suggest that the SZ peak coincides with the peak in the pressure times the line-of-sight characteristic 
length and not the pressure maximum (as it would for clusters in equilibrium). The peak in the X-ray emission, 
as expected, coincides with the density maximum of the main cluster. As a consequence, the morphology of 
the SZ signal and therefore the offset between the SZ and X-ray peaks change with viewing angle. As an ap- 
plication, we compare the morphologies of our simulated images to observed SZ and X-ray images and mass 
surface densities derived from weak lensing observations of the merging galaxy cluster CL0152-1357. We find 
that a large relative velocity of 4800 km/s is necessary to explain these observations. We conclude that an anal- 
ysis of the morphologies of multi-frequency observations of merging clusters can be used to put meaningful 
constraints on the initial parameters of the progenitors. 

Subject headings: galaxies: clusters: general-galaxies: clusters: intracluster medium -X-rays: galaxies: 
clusters-methods: numerical -galaxies: clusters: individual (CL0 152- 1357) 



1. INTRODUCTION 

Our most successful cosmological model, the cold dark 
matter model dominated by the cosmological constant 
(ACDM), predicts that clusters of galaxies form from the 
largest positive matter density fluctuations. The distribution 
and evolution of these rare, very large positive density fluctu- 
ations are extremely sensitive to the underlying cosmological 
model. Therefore clusters have been extensively used to put 
constraints on cosmological models. In general, we can de- 
fine two main categories to determine cosmological parame- 
ters using galaxy clusters: individual and statistical methods. 

Individual methods use accurate measurements of clusters 
to derive cosmological parameters. The most extensively 
used individual method is the SZ-X-ray method, which 
is taking advantage of the different dependence of the 
SZ signal and the X-ray emission on the physical param- 
eters of galaxy clusters to derive the distances to them 
directly. Measuring the redshift to clusters, this method 
is making use of the distance- redshift function t o de- 
rive cosmological parameters dBonamente et alJ 120061 : 
Schmidt et al. 2004; for prospects for future surveys see 
Molnar et al J 12004 iM olnar, Birkin shaw & Mushotzk"vll2002l: 
Haiman. Mohr & Holder! 120011: iHolder. Haiman & MoM 
20011) . Another individual method, for ex ample, is based on 
the gas mass fr action-redshift function dAllen et alJ 120081 : 
lEttori et ail 120091) . In theory, resonant line scattering along 
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with X-ray or SZ imaging can also be used to derive distances 
to clusters (Molnar, Birkinshaw & Mus hotzky|2006l) . 

Statistical methods use the average properties of many 
clusters to derive cosmological param eters, such a s, for 
example, the cluste r mass function dMantz et alJ 120101 : 
Vikhli nin et alJ 120091). the X-ray l umino si ty and temper- 
ature function dDelPopoloetalJ l20lot Schuecker et al. 

lAngri ' 



20031: iBorgani et all 120011: lAngrick & Bartelmannl 
201 lb . Using SZ and X-ray number counts based 
on future surveys are also a promising pos s ibility 
(iMaiumdar & Mohrl 120031: iLevine. Schulz & Whitel 120021: 
i Weller, Battve & Kneissil 120021; iHaiman, Mohr & Holder! 
200 U IHolder, Haiman & Mohrj 12001b . Recent re- 
views on using galaxy clusters for cos molog y can 
be found in [Allen, Evrard & Mantzl (120111) and 
iRosati. Borgani & Normanl (120021) . 

For both main methods of using clusters to determine cos- 
mological parameters, it is crucial to understand the physics of 
clusters, the distribution of the different components, and their 
formation history. Our ACDM models predict that massive 
clusters form by merging. Therefore some fraction of clus- 
ters must be in a merging state. Merging clusters with small 
mass ratios result in enhanced X-ray luminosities and tem- 
peratures relative to relaxed clusters, as shown by using nu- 
merical simulations (Randall et al. 2002; Ritchi e & Thomas! 
120021; iRicker & Sarazinll2001l) . The derived mass for a clus- 
ter in this stage will be overestimated, thus the derived mass 
function will be biased. When using individual methods it is 
crucial to be able to identify mergers and exclude them from 
our analysis. Statistical methods need to be corrected for the 
affect on different stat istical properties o f clusters. Assuming 
a LCDM cosmology, IWik et alJ (12008b have found that the 
maximum of the Comptonization parameter, y, is increased 



2 



Molnar, et al. 



substantially during the first core passage, resulting a signifi- 
cant bias of 20%-4Q% in determining £!,„ and erg, while cos- 
mological parameters derived from the integrated Compton-y 
parameter introduce only less or equal to 2% bias. 

Subarcminute resolution Sunyaev-Zel'dovich (SZ) effect 
images of some clusters of galaxies compared with high reso- 
lution X-ray images show that the positions o f the maxima of 
the SZ and X-ray signals d i ffer significantly dMassardi et al. | 
' 2010t IKorngut et all I I2010t iMason et al.l 120 lOt iMalu et al.1 
2010t IRodriguez-Gonzal vez et al.1l2010t) . The observed off- 
sets imply that these clusters are not in dynamical equilib- 
rium, therefore the derived physical parameters fro them, and 
the derived cosmological parameters would also be biased. 

In this paper we study quantitatively the offsets between 
the SZ and X-ray peaks after the first core passage using 
3-dimensional N-body/hydrodynamical simulations of ideal- 
ized binary merging clusters of galaxies using FLASH. We 
focus on the offsets after first core passage since this phase is 
the easiest to analyze both observationally and theoretically. 

2. OFFSETS BETWEEN SZ AND X-RAY PEAKS IN GALAXY 
CLUSTERS 



Massardi etaT] d2010l) observed CL JO 152- 1347 using the 
Australia Telescope Compact Array ( ATCA) at 1 8 GHz with 
an angular resolution of 35" x 35". CL JO 152- 1347, at a 
reshift of 0. 83, is one of the most massive high redshift clus- 
ters known. Massar di et al.1 (|2010) have found that the X-ray 
center based on XMM-NEWTON and their ATCA observa- 
ti ons are disp l aced b y about 342 kpc («45"). 

IMalu et all (120101) observed the "bullet cluster " (IE 0657- 
56) at a resolution of 30" using ATCA. They found two X-ray 
peaks at about 1.5'apart at the East (E) and West (W) part of 
the cluster. The peak at W belongs to the infalling cluster, the 
"bullet". The offset between the SZ and X-ray peaks were 
found to be about 35". 

An offset of about 20" was found between the S Z and 
X-ray pe aks in RX J 1347- 1145 by IKorngut etaD d2010h 
(see also Mas on et al.l 120 lOl) using Mustang on GBT with 
about 10"— 1 8" resolution confirming an offset which had 
been discovered earlier using the lower resolution and sen- 
sitivity No beyama Bolometer Array on the Nobeyam a 45-m 
telescope (Komatsu et al. 2001; iKitavama et al. 2004). Using 
Mustang, Korng ut et al.l (|2010) observed two more disturbed 
clusters and found a kidney-shaped SZ feature between the 
two peaks of the X-ray emission displaced by about 20" in 
MACS0744, and a highly disturbed SZ distribution with mul- 
tiple peaks with a prominent ridge like feature oriented per- 
pendicular to the line connecting the center of the main and 
secondary total mass distributions in CL1226. Their prelimi- 
nary study of disturbed clusters demonstrated the importance 
of high-resolution SZ observations in identifying shocks in 
g alaxy clusters particularly at high r edshifts . 

IRodriguez-Gonzalvez et al.l (120101) found an about 20" dis- 
placement between the X-ray and SZ peaks in A2146 using 
the Arcminute Microkelvin Imager (AMI). AMI has a 3' res- 
olution in its most compact configuration, which is suitable 
for extended sources, and an about 30" resolution in its more 
extended configuration to remove point sources. 

3. FLASH SIMULATIONS OF MERGING GALAXY CLUSTERS 

A number of self-consistent 3-dimensional (3D) binary 
merger simulations have been carried out recently using La- 
grangian (Ritchie & Thomas 2002; McCarthy et al. 2007; 
Poole et al. 2006, 2007), and Eulerian codes (Ricker & 



Sarazin 2001; ZuHone 2010) to study the effects of differ- 
ent mass ratios, impact parameters and cluster models on X- 
ray emission, SZ signal and mass surface density. Binary 
merger simulations were also carried out to compare numer- 
ical simulations with X-ray and SZ observations of individ- 
ual clusters (for the "bullet cluster", CI 1E0657-56, Springel 
& Farrar 2007 and Mastropietro & Burkert 2008; and for 
Cl 0024+17, ZuHone et al. 2009a,b). Binary merger sim- 
ulations using self consistent N-body smoothed particle hy- 
drody namics (SPH) and Eul erian (FLASH) have been carried 
out by Mitchell et al. (2009|) to compare the r esults for mix- 
ing/tur bulence for these two different methods. M itchell et al.l 
(2009) showed that SPH codes suppress turbulence, while 
adaptive mesh refi nement (AMR) cod es treat them more re- 
alistically (see also Agertz et al. 20071). 

Since turbulence is important in cluster mergers, we choose 
to use a publicly available parallel Eulerian parallel code, 
FLASH, developed at the Center for Astroph ysical Ther- 
monu clear Flashes at the University of Chicago dFrvxell et al.l 
2000). FLASH is using the Piecewise-Parabolic Method 
(PPM) of Colella & Woodward (1984) to solve the equations 
of hydrodynamics, and a particle-mesh method to solve for 
the gravitational forces between particles in the iV-body mod- 
ule. The gravitational potential is calculated using a multi- 
grid solver (Ricker 2008). FLASH uses adaptive mesh re- 
finement (AMR) with a tree-based data structure allowing re- 
cursive grid refinements on a cell-by-cell basis on a Cartesian 
grid. 

Our simulations achieved a 12.7 kpc resolution at the clus- 
ter centers and the merger shocks. Our box size, 13.3 Mpc 
on a side, was large enough, therefore there was no need for 
corrections for mass loss. 

3.1. Initial Conditions 

We assumed spherical cluster models with a cut off of 
the distrib ution of the dark matte r and gas at the virial ra- 
dius. R v ,r dBrvan & Normanll 19981) . We used an NFW model 
dNavarro. F renk & White fl997l) for the dark matter density, 



PDM(r) = 



x(l +x) 2 ' 



(1) 



where x = r/r s , and p s , r s = r V i r /c v i r are scaling parameters for 
the density and radius, c V! > is the concentration parameter, and 
r < Rvu- The gas distribution was assumed to be a truncated 
non-isothermal j3 model, 



p(r) = 



Po 



(1 + 3,2)3/3/2 < 



(2) 



where y = r/r core , and po, r core are the central density and gas 
scale radius, and r < R v i r . The temperature of the gas was 
determined form the equation of hydrostatic equilibrium via 
numerical integration. The equation of state for the gas was an 
ideal gas equation of state with 7 = 5/3, and the mean atomic 
mass was p = 0.592. We adopted r core = 0.12 /? v ; r , and (3=1, 
which are consistent with our analysis of numerical simula- 
tions dMolnar et al.1 120101) . In our simulations presented in 
this paper, we used 5 and 8 for c V! >, for the main and infalling 
(sub)cluster following the trend that less massive clusters are 
more concentrated. We have also run some simulations with 
different concentration parameters to check their effect. As- 
suming M,„, and the gas mass fraction of 0. 14 we derive /?„>, 
Ps, r s , po. We treat the small fraction of baryonic matter in 
galaxies along with the dark matter since they can be assumed 
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TABLE 1 
Initial Parameters 



ID 


Mi 


M 2 


V 


P 


RMlV48p()0 


2.1 


1.0 


4800 





RMlV48plO 


2.1 


1.0 


4800 


100 


RMlV48pl5 


2.1 


1.0 


4800 


150 


RMlV48p20 


2.1 


1.0 


4800 


200 


RMlV48p25 


2.1 


1.0 


4800 


250 


RMlV48p35 


2.1 


1.0 


4800 


350 


RMlV45pl5 


2.1 


1.0 


4800 


350 


RMlV40pl5 


2.1 


1.0 


4000 


150 


RMlV35pl5 


2.1 


1.0 


3500 


150 


RMlV30pl5 


2.1 


1.0 


3000 


150 


RMlp6V48P15 


2.1 


1.6 


4800 


150 


RMlp3V48P15 


2.1 


1.3 


4800 


150 


RM0p7V48P15 


2.1 


0.7 


4800 


150 



NOTE. — See text for other parameters 

to be collisionless for our purposes. Therefore our dark mat- 
ter particles also represent baryonic matter in galaxies. The 
number of dark matter particles at each cell was determined 
by the density and the total particle number (5 million). 

The velocities of the dark matter particles were determined 
by sampling a Maxwellian distribution (local Maxwellian ap- 
proximation) wi th the velocity dispersio n, er r derived from the 
Jeans equation (Lokas & Mamon 2001). Assuming isotropic 
velocity dispersion (the angular and radial components are 
equal: ag = o>), the Jeans equation becomes: 



1 d 
p dr 



d<l> 



(3) 



where <£>(r) is the gravitational potential, which, in our case of 
an NFW distribution, becomes 



V and P represent the mass of the second cluster in units of 
10 14 M Q , the initial relative velocity in units of 100 kms" 1 , 
and the impact parameter in units of 10 kpc. The mass of the 
main cluster was assumed to be 2.1 x 10 14 M Q , The concen- 
tration parameters were assumed to be 5 and 8 for the main 
and the infalling cluster. We have also run a few simulations 
with different gas mass fractions, concentration parameters 
and mass rations, and one simulation with higher resolution. 

4. SZ, X-RAY AND SURFACE MASS DENSITY IMAGES 

After the simulation finished, we generated SZ, X-ray and 
surface mass density, E, images assuming different viewing 
angles (expressed as rotation angles out of the plane of the 
sky) and phases (time elapsed after first core passage) of the 
collisions. 

Since relativistic corrections are important in the high tem- 
perature shocked gas (T s > 30 keV), we generated SZ and X- 



ray images using relativistic corrections. 
SZ surface brightness using 



We calculated the 



Isz(x,y) oc 



/ 2 p g T g [g(u) + E^Y n e"]d£, (7) 
Je, 



where O = kgTg / (m e c 2 ) and g{v) = coth(jc„/2) - 4, x v = 
hpv/QiBTcmb), and v is the frequency, T cm b is the temperature 
of the cosmic microwave background, hp, kg, c and m e , are 
the Planck and Boltzmann constants, the speed of light and 
the electro n mass. The relativistic correctio ns, Fi.2,3.4, were 

taken from lltoh. Kohvama & Nozawal dl998h. 

We generated the X-ray images using lRvbicki & Lightmanl 
(1979)'s expression for the relativistic X-ray thermal 
bremsstrahlung, 



$(r) = -V v 2 g(c vir ) 



Inq 



(4) 



where the circular velocity is V~ = GM vir /R vir , 



g(c vir ) = [ln( 1 + c,„>) - c vir /( 1 + c vir )} 1 , (5) 

and q = 1 +xc vlr . Thus, from Equation[3] we obtain the veloc- 
ity dispersion: 



of(r) = V v 2 g(c vir )xq 2 





ln<? c vir 


f 


X^q x q' 



dx, 



(6) 



(see lLokas & Mamon 2001 for more details on the properties 
of NFW profiles). The direction of the velocities were as- 
sumed to be isotropic. Running a control simulation with an 
isolated cluster, we found that this local Maxwellian approxi- 
mation will relax the density distribution close to the assumed 
NFW model (as expected). Therefore we conclude, that this 
approximation is adequate for our purposes, since we are in- 
terested in the effect of merging on the gas distribution and not 
in the details of the changes in the dark matter distribution. 

3.2. FLASH 3D Simulations 

We have run a set of simulations of galaxy cluster mergers 
systematically changing the initial mass ratios, impact param- 
eters, and relative velocities. In Table Q] we show the initial 
parameters for those runs we discuss in detail in our paper. 
In this table the first column is the identification number for 
our runs using the following convention: the numbers after M, 



h(x,y) oc 



2 T 



1/2 



gff (1+4.4 x l0- l0 T g )d£, (8) 



where the Gaunt factor, g ff = (2 V(3)/7r)[ 1 + 0.79(4.95 x 

10 5 /7p] (see also lHughes & BirkinshawfT^l) . 

We integrated the total (dark matter and gas) density along 
the line of sight (LOS) from different viewing angles to obtain 
the mass surface density images at position, x and y: 



T, x (x,y)oc / (p d + p g )d£, (9) 
where I is the spatial coordinate in the LOS. 

5. RESULTS: OFFSETS BETWEEN THE X-RAY AND SZ CENTERS 

In this paragraph we discuss general aspects of our results. 
We have run simulations with different initial galaxy cluster 
masses and found similar dependence of the offsets between 
X-ray and SZ centers on different impact parameters and rel- 
ative velocities. Thus we choose to show our results with 
fixed Mi = 2.1 x 10 14 M Q , since these give the best match 
with the observed morphology of CL0 152- 1357, which we 
will discuss in the next section. In Figure Q] we show snap- 
shots of merging galaxy clusters after the first core passage 
assuming different impact parameters, P = 0, 150 and 350 kpc, 
all other parameters of the simulations were held fixed (runs 
RM0p7V48P15 and RMlp3V48P15, see Table 1 for details of 
the simulations). The contour plots centered on the main clus- 
ter show the mass surface density, S, the X-ray emission and 
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FIG. 1. — Snapshots of merging galaxy clusters (infall velocity is upward) after the first core passage at time, T, in output time units (left to right), assuming 
different impact parameters, P = 0, 150, and 350 kpc (from top to bottom). All other parameters held fixed (runs RMlV48p00, RMlV48pl5, and RMlV48p35; 
see Table 1 for details). We show contour plots of projected mass surface density, X-ray emission and SZ contours (red, blue and green lines). The contour levels 
are chosen arbitrarily to make it easier to compare morphologies and avoid clutter. The collision is in the plane of the sky (the projection axis is in the line of 
sight). The images are 2.5 Mpc X 2.5 Mpc. 
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FIG. 2. — Same as Figure[T]but assuming different initial relative velocities, 3500 and 4000, in unit of km s 1 (from top to bottom), all other parameters held 
fixed (runs RMlV35pl5 and RMlV40pl5; see Table 1 for details of the simulations). 

while the gas gets shocked because of the collision with the 
gas of the main cluster, and it is slowed down due to ram pres- 
sure. As a consequence, the gas of the infalling cluster will be 
out of equilibrium, and it gets displaced relative to the dark 
matter, it falls behind the center of the infalling cluster dark 
matter. 

The effect of the ram pressure as a function of the impact 
parameter can be seen clearly. At larger impact parameters the 
core of the infalling cluster goes through less dense regions of 
the main cluster and, since the ram pressure is proportional 



SZ contours (red, blue and green lines) projected to the main 
plane of the collision containing the two mass centers and the 
relative velocity vector (i.e., the collision is in the plane of the 
sky). T represents the elapsed time in output file time units. 
The images are 2.5 Mpc x 2.5 Mpc. 

The most obvious feature we recognize looking at Figure Q] 
is the different behavior of the infalling cluster gas and the 
most massive component, the dark matter. As we expected, 
the dark matter in the infalling cluster, due to its collision- 
less nature, simply passes through the main cluster smoothly, 
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T=28 ■■ T = 30 " T=32 " T=34 " T=36 

Ml =2.1 M2 = 0.7 :: M1 =2.1 M2 = 0.7 : M 1 =2.1 W2 = G.7 - M'=2.1 M2 = 0.7 :: M 1 =2.1 M2 = 0.7 




°- b ~ T=28 - : T = 30 :: T=32 T=34 " T=36 ■ 

M 1 =2.1 M2=1.3 :: M 1 =2.1 M2=1.3 : M1 =2.1 M2=1.3 : M1 =2.1 M2=1.3 :: M 1 =2.1 M2=1.3 
o.o r t t i t • 

0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 2.5 

FIG. 3. — Same as FigurefT]but assuming different mass rations: the mass of the main cluster was fixed at Mi = 2.1, and the infalling cluster with masses: M2 
= 0.7 and 1.3 in unit of 10 14 Mq (from top to bottom), all other parameters held fixed (see text for details). 




FIG. 4. — Snap hots of mass surface density, SZ effect and X-ray images of merging clusters of galaxies as a function of different rotation angles, = 0°, 
30° , 60° (bottom to top), away from the main plane of the collision with axis perpendicular to the direction of the initial infalling velocity for different initial 
velocities. The color code is the same as in FigurefT] The infall velocity points upward. The initial masses are: M\ = 2.1 X 1O 14 M0, Mi = 1.0 X 10 14 Mq, 
except for the right column. From left to right: P = 100, 150, 150, 150 kpc, V = 4800,4800,3500,3500 in km s -1 . All other parameters help fixed (see text for 
details). The images are 2.5 Mpc X 2.5 Mpc. 
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FIG. 5. — SZ-X-ray offsets as a function of the distance between the 
dark matter mass centers of the two merging clusters for different initial 
relative velocities: 3000, 3500, 4000, 4500, and 4800 kms" 1 (plus signs, 
stars, squares, triangles and diamonds connected by short dashed, dot-dashed, 
solid, long dashed, dot-dot-dot-dashed lines). 




Moss Peak Distance/Arc Second 

FIG. 7. — Same as Figure[6]but for different initial mass ratios: Mi = 2.1, 
fixed, M9 = 0.7,1.0,1.3,1.6 in units of 1O 14 M0 (stars, squares, triangles and 
diamonds connected by short dashed, dot-dashed, solid, long dashed, dot-dot- 
dot-dashed lines). The plus signs represent results from a higher resolution 
run with Mi =2.1xlO 14 M and M 2 = lxlO 14 M . 
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FIG. 6. — Same as Figure [3] but for different impact parameters: P = 0, 
150, 250, 350 kpc (stars, squares, triangles and diamonds connected by dot- 
dashed, solid, long dashed, dot-dot-dot-dashed lines). All other parameters 
help fixed (see Table 1 for details of the simulations). 

to density x velocity 2 , the infalling cluster retains more and 
more of its gas. The displacement between the SZ and X-ray 
peaks can be seen in all cases with non-zero impact parameter. 

We can see how the displacement depends on the initial rel- 
ative velocity by studying Figure [2] In this figure we show 
results for V = 3500 kms" 1 and 4000 km s" 1 , all other param- 
eters held fixed (runs RMlV35pl5 and RM1 V40pl5; see Ta- 
ble 1 for details of the simulations). Results from simulations 
with the same initial conditions but with relative velocities of 
V = 4800 kms" 1 are shown in the second row of Figure Q] 
Again, we can see that assuming larger relative velocities, the 
infalling cluster can retain less of its gas due to ram pressure 
stripping. Also, larger velocities result larger displacements 
between the SZ and X-ray peaks. However, it seems that if 
the velocity is large enough, ^ 3500 kms" 1 , we obtain large 
displacements, £ 200-300 kpc (the largest offset is at T = 40 
for this relative velocity). 

In Figure [3] we show snapshots of the collisions as a func- 
tion of different masses of the subcluster. The mass of the 
main cluster was fixed at Mi = 2.1 x 10 14 M , and the masses 
of the infalling cluster were M 2 = 0.7 and 1.3 x 10 1 4 M (for 
M2 =1.0 see the second row of Figure Q] for which the other 
initial parameters are the same). Again, all other parameters 
held fixed (see Table 1 for details). In this case we see the 
result of two opposing effects: the ram pressure is striping 
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FIG. 8. — Offsets between the two mass centers and the SZ and X-ray peaks 
(upper and lower set of curves) as a function of the cosine of the rotation 
angle away from the main plane of the collision with axis perpendicular to 
the direction of the initial infalling velocity for different impact parameters: 
100, 150, 200, 250 kpc (stars, squares, triangles and diamonds connected by 
dot-dashed, solid, dashed, dot-dot-dot dashed lines). 

the gas of the infalling cluster, and the gravity is trying to 
keep the gas, the relative velocities were held fixed. Assuming 
M2 = 0.7 x 10 14 M©, the infalling cluster looses all its gas, and 
with increasing mass, the subcluster retains more and more of 
its gas (M 2 = 1, and 1.3 x 10 14 M Q ). 

In Figure |H we show the projections rotated by angles = 
0°, 30°, 60°, away from the main plane of the collision with 
axis perpendicular to the direction of the initial infalling ve- 
locity (from bottom to top, runs RMlV48plO, RM1 V48pl5, 
RMlV35pl5 and RM0p7V48P15). The infall velocity points 
upward. The initial masses are: Mi =2.1 x 10 14 M Q , M 2 = 1.0 
x 10 14 M©, except for the right column for which M 2 = 0.7 
xl0 14 M Q . From left to right: P = 100,150,150,150 kpc, 
V = 4800,4800,3500,3500 in km s" 1 . All other parameters 
help fixed (see Table 1 for details). It is interesting to note the 
significant changes in SZ morphology, and, as a consequence, 
the offset between SZ and X-ray peaks: the offset gets smaller 
with larger rotation angle for all cases shown. 

We show our results for the offsets between the SZ and 
X-ray peaks for different initial velocities, impact parameters, 
and masses in Figure (runs RMlV30pl5, RMlV35pl5, 
RMlV40pl5, RMlV45pl5, and RM1 V48pl5), Fig- 
ure (runs RMlV48p00, RMlV48pl5, RMlV48p25, 
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and RMlV48p35), and Figure Q (runs RM0p7V48P15, 
RM1 V48pl5, RMlp3V48P15, and RMlp6V48P15). 

A qualitative analysis of the offset between the SZ and 
X-ray peaks shows that 40" or larger distances between the 
peaks can be produced as long as the relative velocities are 
at least 4000 kms" 1 (see Figure |5). At large initial relative 
velocities (V = 4800 kms" 1 ) runs with impact parameters be- 
tween 100 and 250 kpc, or, with fixed P = 150 kpc, all masses 
with fixed Mj = 2.1 x 10 14 M Q , and M 2 = 0.7-1. 6 x 10 14 M , 
the offset is larger than 40" (see Figures [6] and |7|. 

We have carried our a run using higher resolution to check 
if our resolution is high enough. The offsets for this run are in- 
cluded in our Figure [7] (plus signs). In this figure, the squares 
represent offsets based on a run we used in all simulations. 
The two offsets are almost identical, thus we conclude that 
our resolution is sufficient for our purpose. 

In Figure [8] we illustrate how the offsets between the SZ 
and X-ray change with rotation angle (out of the plane of the 
collision), 9, for different impact parameters, P = 100, 150, 
200, 250 kpc (runs RM1 V48pl0, RM1 V48pl5, RM1 V48p20 
and RM1 V48p25). The upper lines represent the distances 
between the two mass peaks, the lower lines show the off- 
sets between the SZ and X-ray peaks. Assuming = 0, we 
find large distances between the projected mass centers, about 
1 10", when the largest offset between the SZ and X-ray peaks 
are about 40". At each O, the largest offset we obtain is for 
P= 150 kpc, Mi =2.1 xlO 14 M , M 2 = 1.0 xlO I4 M Q , and 
V = 4800 kms -1 . This is the largest offset for all masses we 
considered here. 

We conclude that large relative velocities of V it 4000 
kms" 1 with different impact parameters close to the core ra- 
dius, and masses can easily produce 40", or larger offset be- 
tween SZ and X-ray peaks. 

6. APPLICATION TO GALAXY CLUSTER CL0152-1357 

In this section we use our results from FLASH simula- 
tions to interpret multi-wavelength observations of the galaxy 
cluster CL0152-1357. We use the morphology of the X-ray, 
SZ and S images, and the offsets between the SZ and X- 
ray peaks as well as the distances between the mass peaks of 
the two components to constrain the initial conditions of the 
collision. We also discuss the implications of our results for 
CL0152-1357 as a test for ACDM and the effect of SZ reso- 
lution on the interpretation of galaxy clusters. 

6.1. Constraining Initial Conditions 

The ATC SZ image of CLO 152- 13 57 seems to be nearl y 
circularly symmetric (see Figure 1 of lMassardi etalj f2010). 
Searching for similar SZ morphology in our projected images 
we find candidates with different impact parameters and ini- 
tial relative velocities (see green contours in the left column 
of Figure HI results for runs RMlV48p25, RMlV35pl5 and 
RM1 V48pl5), also with different initial masses (green con- 
tours in first panel in Figure [9] run RM0p7V48P15 with Mi = 
2.1 x 10 14 M Q and Mi = 0.7 x 10 14 M (see Table 1 for other 
parameters). We assumed a 35" resolution for SZ observa- 
tions (ATCA), and choose the co ntours to reflect the s ignal 
observed using ATCA (Figure 1 of Massardi et al.1120101) . For 
comparison, we also plot the centers of the X-ray emission 
to show that, in all of these cases, there is an offset between 
the SZ and X-ray peaks. These panels illustrate that an ob- 
served nearly circular SZ distribution, with the assumed SZ 
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FIG. 9. — Contours of SZ, X-ray and mass surface density (green, blue, 
red colors) for different models. Left column: similar SZ morphology to 
the observed (to aid the eye, we included the center of the X-ray emission) 
with different initial velocities, 3500, 4800 in units of kms -1 , and impact 
parameters (150 kpc and 250 kpc). Right column: top 2 panels: SZ and X- 
ray morphology similar to the observed, bottom panel: SZ, X-ray and mass 
surface density morphology similar to the observed (60° rotation away form 
the plane of the collision). The contour levels are chosen arbitrarily to make 
an easy comparison to the observed morphology. The images are 350" X 
350". 

sensitivity and resolution, can not be used to identify a shock 
in a cluster, nor can the initial conditions be constrained using 
only SZ observations if we know that there is a shock. It is 
clear that we need a multi-frequency approach in this case. 

Including the observed SZ and X-ray morphology of 
CL0 152-1357 (Figure 1 o f iMassardi et ai1 l20Tol and Figure 
3 of iMaughan et alJl2006h to constrain the initial conditions, 
we find that the best match with observations is provided by 
run RM0p7V48P15 with masses of 2.1 xlO 14 M and 0.7 
x 10 14 M , and 6 = (the collision in the pane of the sky), P 
= 150 kpc, and initial relative velocity 4800 kms" 1 (top panel 
in the right column of Figure [9). We assumed a 12" X-ray 
resolution (XMM-Newton). This high velocity seems to be 
needed to obtain the large offset, about 45" between the SZ 
and X-ray peaks, and the observed SZ and X-ray morphol- 
ogy. Allowing for some uncertainty in the relative pointing 
between the SZ and X-ray instruments, we find a good match 
with observations for run RMlV48pl5 with 2.1 xlO 14 M 
and 1 x 10 14 M Q , and the same impact parameter and relative 
velocity as the previous run (middle panel in the right column 
of Figure |9). This run with a rotation angle of O = 60° can 
also reproduces the observed SZ and X-ray morphology, but 
the SZ and X-ray offset is only about 20". Note, that the rel- 
ative coordinate positions of the SZ and X-ray peaks differ 
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for these two cases, but that can not be observed using these 
two wavelengths because we do not know the original orien- 
tation of the two clusters. Even though the possible initial 
configurations are more restricted by using X-ray and radio 
observations, we can conclude that, using morphology, these 
observations are not sufficient to confine the initial conditions 
of the collision. 

We include now the mass surface density map of CL0152- 
1357 derived from gravitational lensing obs ervations in our 
analysis assuming a resolution of 20" dJee et alJl2005l) . The 
importance of the lensing observations is that the distance be- 
tween the two dark matter centers constrain the phase of the 
collision, i.e. the time elapsed since the first core passage 
(subject to projection effects). Having some constraints on 
the time elapsed since the first core passage and the X-ray 
and SZ morphology enable us to put more constraints on the 
initial conditions of the collision. 

We are looking for an offset between the SZ and X-ray 
peaks of about 45", and a distance of about 50" between the 
mass peaks, in projection. For all projections perpendicu- 
lar to the main plane of the collision (8 = 0"), the distance 
between the mass peaks are around 100"-120", so we need 
to rotate the infalling cluster center towards the line of sight, 
thus we conclude that the rotation angle, 6, should be greater 
than zero. We find that, with 6 = 60°, our run RMlV48pl5 
with P = 150 kpc, Mi = 2.1 x 10 14 M Q , M 2 = lxlO 14 M , 
and V = 4800 kms" 1 , provides the best fit to all three ob- 
servations. The morphology of this model is close to the 
observed SZ and X- ray morphologies (Massardi et al.l l2010t 
Maug han et al.| [2006). and the distance bet ween the projec ted 
mass centers is the same as the observed dJee et al.ll2005l) . It 
is encouraging that these mass es are similar t o the masses de- 
rived for these components by (Jee et al. 2005). We were run- 
ning simulations with larger masses (Mi > 2.1xlO 14 M ), 
but the resulting X-ray morphologies did not match the ob- 
servations. The offset between the SZ and X-ray peaks is 
less than the observed 45". This difference may be caused by 
absolute calibration for the positions between the ACT and 
XMM-NEWTON. However, it is also possible that the lens- 
ing observations have some systematics. Gravitational lens- 
ing observations are very difficult, and the mass surface den- 
sity reconstruction depends strongly on the assumed priors. 
Compare, for example, the results using different priors for 
reconstruction of the mass surface dens ity for RX J 1347.5- 
1 145 (see Figure 4 of Mir anda et al.l2008h . It would be useful 
to have mass surface density reconstructions using different 
methods to check for systematics in the analysis of CL0152- 
1357. 

6.2. Tests for ACDM 

Merging galaxy clusters can provide tests for our most 
successful cosmological scenario, the ACDM models many 
ways. The features of mergers predicted by ACDM models, 
for exa mple, can be compared to those derived from obser- 
vation s dForero-Romero et ail 120101; iMastropietro & Burkerj 

S: tJ -^ 

iForero-Romero et alj d2010l) derived the distribution of the 
offsets between dark matter (DM) and gas centers of galaxy 
clusters drawn from one of the largest non-radiative, ACDM 
cosmological simulations using SPH, usually referred to as 
the MareNostrum Universe. They found that the distribution 
of the 2D offsets between the DM and gas density peaks of 
the bullet cluster would have a probability of l%-2% from 



simulations. Thus they concluded that it is possible to repro- 
duce the observed offset between DM and gas density peaks 
as large as observed in the bullet cluster assuming ACDM 
cosmology. 

The estimated initial relative veloc ities of clusters can also 
be use d as a test for ACDM models. Mastro pietro & Burkertl 
(2008) carried out galaxy cluster merger simulations and 
found that their best model requires a relative initial veloc- 
ity of 3000 kms" 1 for the "bullet" cluster, 1E0657-56. How- 
ever, they concluded that the probability of finding such a 
high relative velocity given their best model is less then 0.5%. 
We found a relative velocity of 4800 kms" 1 for our best 
model for CL0 152- 1357. The median relative velocity of 
a galaxy cluster with about the total mass of CL0 152- 1357 
Pee et all 120051). M i = 5 xl0 14 M Q , is about 1500 kms -1 
faavashi & White! 12006b . The infall velocity of our best 
model is about three t imes the value of our be st fit model. 
Based on the results of Havashi & White (2006), this relative 
velocity is very unlikely in a ACDM model. These results 
suggest that there is a mismatch between the initial relative 
velocities derived from simple idealized cluster merging sim- 
ulations and those derived from ACDM simulations. There 
are two possibilities: either our simple cluster models or our 
ACDM models are missing some important physics. A de- 
tailed analysis of the origin of this mismatch is out of the 
scope of our paper, we are going to address this question in 
a future paper. 

6.3. The Effect of SZ Resolution 

When using clusters to derive cosmological parameters we 
assume that they are relaxed. The statistical methods use 
scaling relations derived based on relaxed clusters, individ- 
ual methods assume relaxed clusters to derive their physical 
parameters. Both statistical and individual methods may be 
biased due to the fact that some clusters are not relaxed. Sta- 
tistical methods might have to be corrected for this effect, 
when using individual methods, we can just simply exclude 
non-relaxed clusters. Some of these aspects of SZ resolu- 
tion and the imp ortance of high-resolution SZ observation s 
were discussed inlKorngut et all d2010l) : lMason et al.l (1 20101) : 
Massardi et al. (2010) for example. 

In this paper we discuss some aspects of the identification 
of non-relaxed clusters focusing on SZ observations. The ef- 
fect of non-relaxed clusters to determining cosmol ogical pa- 
ramete rs using SZ observations was discussed in iWik et all 
(2008). The new generation of ground based SZ telescopes 
have a very large span in resolution, from about 10" to 3'. 

We illustrate the importance of high resolution SZ obser- 
vation in identifying merging clusters in Figure [10] As an 
example, we use the best model we found for CL0152-1357. 
In this Figure we show snapshots for this model after the first 
core passage using different SZ resolution (green contours). 
We also plot the centers of projected mass (red contours) and 
X-ray peak (blue contours). From this Figure we can see that, 
with high resolution, 10"— 35", SZ observations would be suf- 
ficient to identify the merger based on its disturbed morphol- 
ogy. Having also X-ray observations would help us because 
we can use the offset between the SZ and X-ray peaks. In the 
case when the collision is in the plane of the sky (9 = 0), this 
offset is about 50" with high SZ resolution, and about 40" for 
low, 1'- 2' resolution. In this case the offset is observable as- 
suming a better than 40" absolute pointing calibration for the 
SZ telescope, which is a reasonable assumption. However, if 
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FIG. 10. — Projected mass surface density, X-ray emission and SZ contours red, blue and green contours) of run RMlV48pl5 (see TablefT} as a function of 
rotation angle, 0, away from the plane of the collision for different resolutions, FWHM, of the SZ observations: FWHM = 10", 35", 60", and 120" (left to 
right). The rotation angles are: = 0°, 30° , 60° away from main plane of the collision (from bottom to top). The images are 350 " X 350 " . 



the rotation angle is large, say 60°, high resolution SZ instru- 
ments still measure an about 40" offset between the SZ and 
X-ray peaks, but with lower resolution, we obtain an only 10" 
offset, which is close or less then the absolute pointing cali- 
bration of SZ instruments. Therefore, we conclude that with 
an SZ resolution of about 10" and 35", GBT/Mustang and 
ATCA have the sufficient resolution to identify this merger, 
but using lower resolution instruments of l'or 2', BOLOCAM 
and AMiBA, will not be able to identify the merger. We could 
identify the merger in any case if we have not only SZ and X- 
ray measurements, but also surface mass density maps from 
lensing optical/infrared (optical/IR) observations. 

Our results suggest that if we want to identify mergers, we 
should use either high resolution SZ observations, or multi- 
frequency observations. In any case, galaxy cluster observa- 
tions with low resolution SZ instruments with large field of 
view are still useful because they can be potentially used to 
constr ain the large scale distribution of the intra-cluster gas 
(e.g. JMomaretaLll20Toh . 

7. CONCLUSION 

We have been studying galaxy cluster merging using self- 
consistent 3D N-body/hydrodynamical simulations of ideal- 
ized clusters. We have found that significant offsets results 
between the X-ray and SZ peaks (>40") if the relative ve- 
locities are equal or larger than about 4000 kms" 1 . This offset 
can be observed using an SZ instrument with high sensitivity 



and high angular resolution (at least about 10"— 35"), in which 
case the merging cluster can be identified using SZ observa- 
tions only (Figure [Toll. However, using multi-frequency ob- 
servations seems to be a more practical way to identify merg- 
ing galaxy clusters. 

We have applied numerical simulations to interpret multi- 
frequency observations, radio (SZ), X-ray and optical/IR 
(lensing reconstruction) of CL0152-1357. We have found that 
the offsets between the SZ and X-ray peaks provide important 
constraints on the initial parameters of the merging clusters, 
the offsets between the two peaks of the mass surface densi- 
ties provide constraints on the phase of the collision (i.e., the 
time elapsed after the first core passage), and that the SZ peak 
coincides with the peak in the pressure times the characteris- 
tic length in the line of sight and not the pressure maximum 
(as it does for clusters in equilibrium). The peak in the X-ray 
emission, as expected, coincides with the density maximum 
of the main cluster. As a consequence, the morphology of the 
SZ signal, and therefore the offset between the SZ and X-ray 
peaks, change with viewing angle. 

We conclude that analyzing the morphology of SZ, X-ray, 
and surface mass density images based on multi-frequency 
observations (radio, X-ray and optical/IR) enables us to put 
meaningful constraints on the masses of the colliding clus- 
ters, the impact parameter and initial relative velocity of the 
collision. 
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